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Abstract. If the cluster formation rate is constant and the disruption time of clusters 
depends on their initial mass as tdis = tiX the values of t4 and 7 can be 

derived in a very simple way from the age and mass histograms of large homogeneous 
samples of clusters with reliable age and/or mass determinations. We demonstrate the 
method and derive the values of ti and 7 from cluster samples in selected fields of 
four galaxies: M51, M33, SMC and the solar neighbourhood. The values of 7 are the 
same within their uncertainties in the four galaxies with < 7 >= 0.62 ±0.06. However, 
the disruption time of clusters of 10^ Mq is very different in the different locations. 
It is shortest in the inner region of M51 (t 4 = 4 x 10^ yrs) and highest in the SMC 
(t 4 = 4 X 10® yrs). 


1 Introduction 

Studies of the cluster formation history in different galaxies, based on observed 
cluster samples, rarely take into account the fact that clusters may have been 
disrupted. We have started a series of studies to disentagle the effects of cluster 
formation and cluster disruption. Here we present the method for the determina¬ 
tion of the disruption time in galaxies with a nearly constant cluster formation 
rate, with an application to selected regions in four galaxies. 

Oort [12] noticed the lack of old open clusters in the solar neighbourhood and 
derived statistically that Galactic clusters disrupt on a timescale of 5 x 10® years. 
Several authors have found that the mean age of clusters in the LMC and SMC 
is larger than that of clusters in the solar neighbourhood and cloncluded that 
the disruption time must be longer in the MCs (e.g. [9]). Actually, the disruption 
time of clusters is expected to depend on their mass and on the local conditions in 
their host galaxy. Modern observations of large homogeneous cluster samples, e.g. 
with HST, enable the determination of the cluster disruption time in different 
galaxies. 

The disruption time of clusters is expected to scale with the relaxation time 

idis ~ trxt (1) 

where Rd is the cluster radius and Md is the cluster mass (see e.g. [13],[14]). 
If all clusters in a limited region of a galaxy had the same radius, then Eq. 
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( 1 ) would imply that tdis ~ ■ If the cluster radius depends on its mass as 

Rc\ = Cgai X then we expect tdis ~ with 7 = (1 + 3(5)/2, where Cgai 

is a constant that depends on the ambient pressure in the region of the galaxy 
where the clusters are located. (Obviously, this expected relation does not hold 
for clusters in highly elliptical orbits.) 

In this paper we derive the empirical relation between the disruption time 
and the mass of clusters in different locations of galaxies, based on the mass 
and/or age distributions of magnitude limited cluster samples. 

2 The predicted mass and age distributions of surviving 
clusters above a given magnitude limit. 

Suppose that: 

1 . the disruption time of clusters depends on their initial mass Md as 

tdis = b X (Md/lO^ (2) 

where td is the disruption time of a cluster with an initial mass of 10 "^ Mq, 

2 . the formation rate of clusters is constant, 

3. the initial mass function of clusters can be written as N{Mc\) ~ ^ with 

a ~ 2 [16], 

4. clusters fade below the detection limit, in the wavelength band that deter¬ 
mines the magnitude limit of the cluster sample, due to the evolution of 
their stars as Fx ~ t~^. The value of C can be derived from cluster evolution 
models, e.g. the Starburst99 models of Leitherer et al. [11]. 

Then it is easy to show that both the age and mass distributions of surviv¬ 
ing clusters above a certain magnitude limit will consist of double power laws 
(shown in Fig. 1) of the type: 


• dN,i/dM,i ~ M/i-“ 

• dN,i/dt - 

• dN^i/dt ~ 


for low mass clusters due to fading 
for high mass clusters due to disruption 
for young clusters due to fading 
for old clusters due to disruption 


These relations are derived by Boutloukos & Lamers [4] (hereafter BL02). 

The vertical shift of the fading lines in Fig. 1 depends on the detection limit 
of the observations and on the number of clusters formed per unit time in the 
region of the host galaxy that is studied. The horizontal shift of the disruption 
line depends on ti. The slope of the disruption line depends on the values of 7 
and a. So both ti and 7 can be derived from the cluster statistics in a simple 
straightforward way. 

These double power laws are strictly valid if the clusters disrupt instanta¬ 
neously when they reach an age of tdis- However, it can be shown that gradual 
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cluster disruption results in similar power law distributions (BL02). In that case 
the transition between the two powerlaws is not a sharp bend, but it is more 
gradual. This is demonstrated in Fig. 2. that shows the predicted mass and age 
distributions in case the cluster mass decreases gradually on a timescale that 
depends on the mass as in Eq. (2). The full lines are the predicted relations for 
gradual disruption, whereas the dashed lines are the relations in case of instan¬ 
taneous disruption. They show a similar double power law as in Fig. 1, but with 
a gradual transition. 



log Mji/Hs 


log t (yrs) 


Fig. 1. Schematic pictures of the predicted mass distribution (left) and age distribution 
(right) of clusters due to fading (full lines) and disruption (dotted lines). The y-axes 
have arbitrary units, (figure from BL02) 




Fig. 2. Predictions for gradual decay of clusters with an age dependent detection limit 
and with a disruption time that depends on the actual cluster mass as given by Eq. (2). 
Left: the predicted mass distribution (full line). Right: the predicted age distribution 
(full line). The dash-dotted lines show the prediction if there was no disruption but 
only evolutionary fading. The dotted lines show the prediction for evolutionary fading 
and instantaneous disruption. The parameters differ from those used in Fig. 1. (figure 
from BL02) 
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3 Application to regions in four galaxies 

We have applied this method to samples of clusters observed with HST in the 
inner region of M51 by Bik et al. [3], in M33 by Chandar et al. [5],[ 6 ],[7] and to 
the cluster samples of the SMC by Hodge [10] and of the solar neighbourhood by 
Wielen [15]. The distributions indeed show the power-law decrease at high mass 
and/or high age due to disruption. Figures 3a shows the mass distribution of the 
M51 cluster sample and Fig. 3b shows the age distribution of the M33 cluster 
sample. (BL02 also studied the age distribution of clusters in M51 and the mass 
distribution of clusters in M33.) The full lines are the predicted distributions for 
evolutionary fading below the detection limit for these samples. (The wiggles at 
low age or small mass are due to the appearance and disappearance of red super¬ 
giants.) At high age or high mass the observed distributions are powerlaws that 
clearly decrease much steeper than predicted for evolutionary fading only. This 
is due to disruption. Figures 4a and 4b show the age distributions of clusters 
in the SMC and the solar neighbourhood repectively. The age distribution of 
the clusters in the SMC and the solar neighbourhood show a flat part at young 
ages, because the clusters do not fade below the detection limit in the samples 
used. Therefore, the distributions at young age show the constant cluster forma¬ 
tion rate. ^From the observed powerlaw relations we derived the values of the 
disruption parameters ^4 and 7 . 

The results are summarized in Table 1. The values of 7 are very similar 
(within the uncertainties) for the cluster samples in the four galaxies, with a 
mean value of < 7 > = 0.62 ± 0.06. The values of ^ 4 , however, are very 

different for the different galaxies: the inner region of M51 has the shortest cluster 
disruption time and the SMC has the longest disruption times. The difference 
amounts to about a factor 10 ^. 


Table 1. The parameters of the disruption time: tais = x (Md/lO"^)^ 


Galaxy 

Region 

Nr 

Age 

Mass 

log ti 

7 


range 

rgai(kpc) clusters log (yrs) 

range 
log Mq 

log (yrs) 


M51 

0.8 - 3.1 

512 

6.0 - 9.7 3.0 - 5.2 7.64 ±0.22 

0.57 ±0.10 

M33 

0.8 - 5.0 

49 

6.5 - 10.0 3.6 - 5.6 8.12 ±0.30 

0.72 ±0.12 

MW 

7.5 - 9.5 

72 

7.2 - 10.0 

— 

9.0 ±0.3 

0.60 ±0.12 

SMG 

0 - 2 

314 

7.6 - 10.0 

— 

9.6 ±0.3 

0.61 ±0.08 






Mean 

0.62 ±0.06 
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Fig. 3. The mass distribution of M51 clusters in the inner 0.8 - 3.1 kpc (left) and the 
age distribution of M33 clusters at 0.8 < roai < 5.0 kpc (right). The full lines are the 
predicted decrease due to fading below the detection limit (from Starburst99 models). 
The dashed lines are powerlaw fits for disruption, (figure from BL02) 



SMC J 


314 clusters : 


Hodge (1987) J 



: —*— 


: — ^ 


: — 


: log dN^bs/dt = +0.68 - 1.648-log(t/10®) 
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Fig. 4. The age distributions of clusters in the SMC (left) and in the solar neighbour¬ 
hood (right). Evolutionary fading below the detection limit is not important for these 
cluster samples, so the young clusters show a constant cluster formation rate. The 
dashed lines are powerlaw fits for disruption, (figure from BL02) 


4 Future studies 

The full study has been described by BL02. Future planned or ongoing studies 
include: 

• study of the disruption time of clusters in different regions of M51 (Bastian 
& Tamers, [1]) 

• study of the disruption time of clusters in the starburst galaxy M82 (de Grijs 
et al. [8]) 

• study of the disruption time of clusters in M81 as a function of mass and 
radius (Bastian et ah, [2]) 

• comparison of the derived disruption times with N-body simulations (Tamers 
& Portegies Zwart, in preparation) 

• comparison of the disruption times of interacting and non-interacting galax¬ 
ies (Bastian et ah, in preparation) 
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